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Abstract. - Motivated by the recent experiment on the Bose-Einstein condensation of 52 Cr 
atoms with long-range dipolar interactions (Werner J. et al., Phys. Rev. Lett., 94 (2005) 
183201), we consider a system of bosons with repulsive nearest and next-nearest neighbor inter- 
actions in an optical lattice. The ground state phase diagram, calculated using the Gutzwiller 
ansatz, shows, apart from the superfluid (SF) and the Mott insulator (MI), two modulated 
phases, i.e., the charge density wave (CDW) and the supersolid (SS). Excitation spectra are 
also calculated which show a gap in the insulators, gapless, phonon mode in the superfluid and 
the supersolid, and a mode softening of superfluid excitations in the vicinity of the modulated 
phases. We discuss the possibility of observing these phases in cold dipolar atoms and propose 
experiments to detect them. 



The recent observation [2] of quantum phase transition of interacting bosons in an op- 
tical lattice has given a new impetus to the study of correlated bosons. The possibility of 
manipulating both the interactions among the constituents as well as disorder in the system 
with extreme control, make them an ideal object over conventional solid state systems to 
study many-body effects. Recent experimental success in obtaining a dipolar condensate of 
52 Cr atoms [1] opens up new directions in the study of lattice bosons with long-range in- 
teractions [3]. This also leads to the possibility of achieving new phases such as insulating 
charge density wave and the supersolid that has both superfluid and crystalline properties [4] . 
The recent experimental observation of nonclassical moment of inertia of solid 4 He indicates 
a possible signature of SS [5], although the existence of the SS is still a controversial issue, 
and we believe, the dipolar condensate in optical lattice can provide a testing ground for it. 

A collection of interacting bosons on a lattice is described by the Bose-Hubbard model [6] 
which captures the physics arising from the competition between kinetic energy of bosons and 
their on-site interaction. However, long-range interactions among the dipolar bosons at differ- 
ent lattice sites make it necessary to include nearest-neighbor interaction in the Hamiltonian. 
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The resulting extended Bose-Hubbard model (eBHM) Hamitonian reads: 

H = afa i+ s + ^^2ni(rii - I) + V^2 ^ih i+ s- (1) 

i, 6 i i,5 

Here creates a boson at site i, hi = af&i is the boson number operator, t, U, V, are the 
nearest-neighbor hopping, the on-site interaction and the nearest-neighbor interaction (NNI) 
respectively, and S represents the nearest neighbors of site i. This model in a high density limit, 
where it corresponds to a quantum phase model (QPM) , has also been studied in the context of 
a wide variety of other problems such as Josephson junction arrays, granular superconductors, 
liquid 4 Hc in Vycor etc. [7,8]. For small t, the role of NNI is to stabilize an insulating state, 
similar to the CDW phase in solids, where the number density alternates at every lattice site 
between two successive integers. In the case of soft-core bosons, it also brings in a new phase 
where both diagonal crystalline order and off-diagonal superfluid long-range order coexist and 
hence called the supersolid [4]; in this phase the number density shows a modulation, while 
still retaining a nonvanishing SF fraction. The model in Eq. |QJ without the NNI has been 
extensively studied previously using different methods [6,9, 10], and the mean-field results were 
found to be in good agreement with the quantum Monte-Carlo (QMC) results [9], though the 
critical behavior is beyond the description of the former. However, one can expect that with 
the increasing dimensionality d of the system, larger filling factors, and the longer range of 
the interactions, mean-field results will progressively become accurate. In this letter we study 
the eBHM using Gutzwiller mean field theory [11] to obtain the ground state phase diagram. 
More importantly, we extend this approach, using the time dependent variational principle, 
to extract the excitation spectrum of the system. This method has the added advantage 
of generalizing to the nonuniform systems, for example in a confining trap, which is always 
present in experiments. We discuss the possibility of observing these novel phases in cold 
atoms and propose ways to detect them. 

In order to find the ground state of the Hamiltonian , we use the variational minimiza- 
tion of H — fj,N with a Gutzwiller wavefunction = Yii J2 n fn" 1 \ni) where \rii) is the 
Fock state with n particles at site i, fi is the chemical potential, and N — J^. hi is the total 
particle number operator. Minimization with respect to the variational parameters fn gives 
the ground state, which could be a MI, SF, SS, or a CDW. First we consider the case when 
Vd < U/2. In Fig. 1, we present the phase diagram for Vd/U = 0.4 [12]. For small enough 
t, incompressible MI and CDW phases are formed. The MI phase with integer no number of 
particles per site (denoted by MI(no)) is characterized by /„ = S n ^ no and a vanishing con- 
densate order parameter 4>i = J^s W ®i+s W- The CDW phase has a two sublattice periodic 
structure such as no, no — 1 particles at the consecutive sites (denoted by CDW(no/2)) and 
zero 4n. The SF and the SS phases have non- vanishing </>; in the former it is uniform and in the 
latter it has a two-sublattice modulation. At small values of /i and t — 0, we have a CDW(l/2) 
with |1, 0, 1, 0) particle density modulation. Upon further increase of /z, MI(1) phase becomes 
energetically favourable. In general, at t = 0, we have transitions between CDW(no/2) to 
MI(n ) at /i = U(n - 1) + 2Vdn and then to CDW(n /2 + 1) at /J, = Un + 2Vdn . As we 
lower the value of V, the area of the CDW lobe reduces continuously and so is the transition 
point at t = between the CDW and the MI phases; the CDW lobe vanishes completely for 
V — 0, thereby reducing to the phase diagram as in Ref. [6]. Uniform, homogeneous SF phase 
appears for sufficiently large values of t. In between the CDW and the SF, a supersolid (SS) 
phase is formed, which has a density modulation as well as nonvanishing <f). Here we must 
stress that a finite onsite interaction is needed, in addition to the NNI, for the existence of 
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checkerboard SS phase which is unstable in the hard-core limit of eBHM [13]. As t increases, 
the system undergoes a second-order transition from CDW(l/2) to a SS for larger values of \i 
(mostly for particle-doping away from half-filling) and a first order transition to SF at smaller 
values of fi (for hole-doping away from half-filling) . A multicritical point exists at fi/U = 0.132 
and td/U — 0.137 where the two second order transition lines between SS-SF and CDW-SS 
meet the first order line between CDW-SF (see Fig. 1). In Fig. 2a, we plot the contours of 
fixed particle density as a function of p and t for Vd/U — 0.4. We find the insulating lobes of 
fixed average commensurate densities no — 0.5, and Uq = 1 etc., and SF and SS phases with 
incommensurate densities (in general). As shown in Fig. 2b and Fig. 2c, we observe at a given 
/x, the density and <fr change continuously across CDW-SS and SS-SF transitions indicating 
a continuous second order transition, whereas a jump in the above quantities signals a first 
order transition between CDW-SF. 

When Vd > U/2, the MI lobes vanish and the new CDW phases such as CDW(l) with 
a particle density modulation |. . . , 2, 0, . . .) appear [14]. At t = 0, the CDW(l/2) transforms 
into CDW(l) when \i > U, which in turn transforms into CDW(3/2) when fi > 2U. In 
general a | . . . , no, 0, . . . > CDW state is formed within a region (no — 1)U < n < noli. With 
increasing t, the supersolid phase is formed around the CDW lobes and at still larger values 
of t, a homogeneous superfluid phase appears. 

In ID due to enhanced quantum fluctuations, our results differ significantly from those 
of the density matrix renormalization group where no SS phase was found [15]. In fact, the 
numerical studies of ID BHM show that the nature of the Mott-SF phase boundary itself 
is quite different from that obtained from the Gutzwiller method. In 2D previous QMC 
studies for the hard-core (U = oo) bosons found unstable checkerboard SS phase which phase 
separates between CDW and SF [13]. Preliminary results for the soft core bosons indicate 
the region of SS phase with negative compressibility [16]. However, a previous QMC study 
of the related quantum phase model by van Otterlo et al. [8] came to the conclusion that a 
finite U (soft core) stabilizes the checkerboard SS phase. Within the Gutzwiller approach, we 
find a stable SS phase over a large region of particle doping, and a relatively smaller region 
of SS upon hole doping the CDW, and a first order transition between the CDW and the 
SF below the multicritical point (region of hole doping) indicating possible phase separation. 
A recent QMC study indeed shows that a stable SS phase exists for soft-core bosons with 
sufficiently large values of V/U and upon particle-doping away from half-filling [17]; however, 
the SS obtained by hole doping turns out to be unstable with respect to phase separation. 
The CDW-SF transition is indeed found to be first order as expected. The region of SS phase 
reduces significantly with reducing V, and we expect that it will be destroyed at small enough 
V due to quantum fluctuations. Increasing the dimensionality to three and/or increasing the 
filling factor will reduce the fluctuations, thus stabilizing the SS further. 

The excitation spectrum of the system can be obtained from the dynamical Gutzwiller 
approach with the variational parameters fn being time dependent. Minimization of the 
effective action zj| — H + fiN \ ^) gives the equations of motion for /„ : 

/W-t^V^+l/^i + ^V^/f i) (2) 

where = J^s n n \fn + ^ I and r is time. The small amplitude fluctuations 5f„\t) around 
the ground state ffi give the excitation spectrum. In the MI ground state, 8fn are connected 
only to <5/,i±i m the Eq. (J2J. The low-lying excitations are creating particles (p) and holes (h) 
in the ground state with a particle density no which needs a finite energy and the system is 



U 



dr 



n (n — 1) — /in + Vnp^ 
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Fig. 1 Fig. 2 



Fig. 1 - Phase diagram of the eBHM for Vd = 0.4 (energies are in units of U). Open circles denote 
the numerical results and the solid lines represent the analytical results. 

Fig. 2 - a) Contour plot of fixed particle density for Vd — 0.4 (energies are in units of U). b) Variation 
of the average particle density for a fixed fi and for Vd = 0.4. Solid line: CDW to SS to SF transition 
at n — 0.3. Dashed line: CDW to SF transition at n = 0.1. c) Variation of 4> for the same parameters. 



gapped. At t = 0, these excitations have energies Uno + 2dVno — n and /i — U(uq — 1) — 2dVno 
respectively. For nonzero t, the dispersions of these excitations are: 

uj ph = ±[-e/2 + U(n - 1/2) + 2dVn - fi] + [(e/2) 2 - eU{n + 1/2) + [/ 2 /4] 1/2 (3) 

where e(k) = 2t Y2t=i d cos (ki)- F° r V = this result agrees with the excitation spectra of the 
MI phase of BHM obtained within the slave boson approach [18]. In general, the CDW state 
has four low- lying excitations in the reduced Brillouin zone (BZ), corresponding to particle 
and hole excitations in each of the A (with n\ particles per site) and B (with ni particles 
per site) sublattices. At t = 0, the energies of particle (hole) excitations of sublattice A are 
E P A = Uni + 2dVri2 — \i, E A — —U(ri\ — 1) — 2dVri2 + [i ,and those of sublattice B are 
E V B = Un,2 + 2dVri\ — fi, and Eg — —XJ(ni — 1) — 2dVn\ + n respectively which, as in the case 
of MI, pick up a dispersion at nonzero values of t. The number of particles on sublattice B is 
ri2 = n% — 1, for Vd < U/2 and n-i = 0, for Vd > U /2. Clearly, the hole excitation in sublattice 
B is unphysical for 112 = 0. For Vd > U/2, there is a qualitative change in the excitation 
spectrum. The excitation corresponding to a particle addition on A becomes energetically 
favourable compared to a particle addition on B. At finite t the excitation spectrum w(k) of 
the insulating states can be obtained analytically by solving the algebraic equation: 

E p A E h A E B E h B - e 2 (k)[K + l)E h A + ^[(na + l)E h B + n 2 E B ] = 0, (4) 

where E A ^ — E P A ^ B ^ — uj. Three low-lying excitations for the CDW(l/2) are shown in Fig. 
3a. The excitations of the MI state can also be obtained from the Eq. © by substituting 
m = ri2 ■ At the phase boundaries of the insulating states excitations become gapless and the 
analytical form of the phase boundaries can also be obtained from Eq. (J3J by substituting 
u) = and k = . These are represented by solid lines along with the numerically obtained 
results in Fig. 1. 

In the SF phase, with the ground state values /A , one can also get the spectra of the 
excitations from Eq. J3J). These quasiparticles, with mixed particle- hole character, are gapless 
and have phonon-like linear dispersion (Fig. 3b). Deep in the SF (for large values of t), 

however, f$ obey the Poissonian distribution fn — <f>f exp ( — \(j)i\ 2 /2 J /\fn\ and Eq. @ 
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Fig. 3 Fig. 4 

Fig. 3 - Excitation spectra for Vd = 0.4 along k x /ir = k y /n = k. a) Excitation spectra of CDW 
(td = 0.08, p = 0.3), solid line: hole excitation on sublattice A, dashed line: particle excitation on 
sublattice B, short-dashed line: particle excitation on A b) Spectrum in the deep SF regime (solid 
line with td = 2, fi = 0.4), dashed line is the excitation spectrum of the effective DNLS for the same 
parameters, c) solid line: the SF excitation spectrum near SS boundary (td = 0.17, pi = 0.4), dashed 
line: the SF spectrum close to CDW phase (t = 0.12, [i = 0.05), notice the mode softening near the 
BZ boundary, d) SS excitation spectrum for td = 0.16, n — 0.4 (energies are in units of U). 

Fig. 4 - Phase diagram of the eBHM in 2D with V = 0.4 and V' = 0.4/(2\/2) (energies are in units 
of U). 

reduces to a discrete nonlinear Schrodinger equation (DNLS): 




In this limit the excitation spectrum can be obtained analytically (dashed line in Fig. 3b) 
with the sound velocity [2td(U + 2Vd)\4>\ 2 ] 1 / 2 and it matches very well with the numerical 
result, obtained from Eq. (J2J). However, as one moves towards the insulating phases, \cj>\ 2 
deviates significantly from p, and the description in terms of the DNLS becomes progressively 
inaccurate. As we approach the MI boundary, <p reduces and finally vanishes at the transition 
point signaling a continuous transition at which a gap opens up for the particle and the hole 
excitations. This behavior has to be contrasted with the approach towards the modulated 
phases (SF-CDW and SF-SS) where a mode softening in the SF spectrum at k = (ir, n) 
implies the instability of the homogeneous SF phase and appearance of states with broken 
translational symmetry [19] (Fig. 3c). Once the SS boundary is crossed, we observe a gapless, 
phonon-like mode in the reduced BZ due to nonvanishing SF fraction, as well as higher energy 
gapped mode as shown in Fig. 3d. The appearance of the gap can be understood as due to 
the doubling of the unit cell and the Bragg reflection at k = (n/2, n/2). 

Let us now discuss the possibility of observing these phases experimentally. The novel, 
translational symmetry broken phases such as SS and CDW appear because of the longer 
range interactions between bosons. There are many ways which we could use to achieve 
bosonic systems with such interactions. One possibility is to use Rydberg atoms; however, 
they have short life-times that make them difficult candidates to observe these phenomena. 
The most promising candidate is a dipolar condensate of 52 Cr atoms which have a large 
magnetic dipole moment [1]. A spin-polarized dipolar system, to the leading order, has a 
power-law tail with Vij ~ ^d/\ r i — r j\ where [id is the dipole moment and i,j are two lattice 
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sites. The on-site interaction U can be tuned by Feshbach resonance and t by the intensity of 
the laser beams which form the optical lattice. Polarization of the dipoles can be manipulated 
by external magnetic field which can produce repulsive interaction between the atoms. Once 
such a condensate is formed in an optical lattice and interactions are tuned by external sources, 
the phases such as the CDW and the SS could be observed in experiments similar to Ref. [2]. 
Both CDW and SS phases are characterized by modulation of the particle density in real 
space imaging (since both of them have a crystalline order). In addition the CDW will have a 
gap in the excitation spectrum. The SS phase will show phase coherence in the time of flight 
experiments. In particular, we find that momentum distribution n(k) in a TV x TV lattice sites 
is given by: 

n(k) = (p A + p B )/2+(\<f) A \ 2 + \<f) B \ 2 )(a x a y -l)/2 + a x a y (\(f) A \ 2 cos(k x + ky) 

+ |0 B | 2 cos(k x - k y ))/2 + 4>A4>BOi x a v {cosk x + cosk y ) (6) 

where a x / y = (2/N)sm 2 (k x / y N/2)/ sin 2 (k x / y ). In the SF, p A — Pb and <f> A = 4>b- We see 
that contrary to the SF case, where the momentum distribution has a sharp central peak, the 
SS is characterized by additional peaks at (±7r, ±7r), (small compared to the central peak) 
characteristic of the ordering wave-vector of the checkerboard phase. We have checked that 
the presence of a confining harmonic trap does not destroy the CDW and the SS phases, but 
leads to coexistence of different phases [20] . Time dependent Gutzwiller theory can be applied 
to calculate the excitation spectra in trapped systems as well [21]. 

Finally, we consider the effect of next-nearest neighbor interaction (NNNI) , 
^' Si S' n i n i+s' ' 011 the phase diagram, where 5 is next-nearest neighbor of site i. In 2D 
for V = Vy(2\/2) (relevant for dipolar interaction), we plot the phase diagram in Fig. 4. 
Checkerboard CDW and SS are formed for V > 2V . In addition new phases appear, such as: 
CDW with one particle (CDW(l/4)) and three particles (CDW(3/4)) in 2 x 2 sublattice and 
SS with three sublattice modulation (SS2) (see Fig. 4). Similar phases can also be obtained by 
mapping hard-core bosons to Hcisenbcrg spin model with NNNI [16,22], although the phase 
diagram is different from the soft-core boson model. In real dipolar gases, the long-range 
nature of the interaction (which is truncated to NNNI in this work) as well as anisotropy 
(which is also neglected here) could be important. However, as can be seen from Fig. 4., 
inclusion of further interactions introduces CDW states of different filling fractions and CDW 
and SS of different ordering wavevectors; we do not expect fundamentally new phases to 
appear due to the long-range interaction. It is also clear that striped CDW or SS phases 
with (it, 0) ordering will be absent in real dipolar gases since 2V < V. Even for realistic 
dipolar interactions, with the anisotropy (e.g. an attractive potential along the z-axis when 
the dipoles are polarized in the x — y plane) included, we expect, for example, CDW and 
SS phases with (ir, n, 0) wavevector, if a collapsed phase is avoided by having a large on-site 
repulsion that can be achieved using a Feshbach resonance or alternatively, by having pancake- 
shaped condensates. This, however, needs further study. Preliminary results indicate that this 
indeed is happening. 

In conclusion, we have studied the lattice bosons with NNI and presented their ground 
state phase diagram and the low-lying excitations in the four phases, by a time-dependent 
generalization of the Gutzwiller ansatz. In particular, the SF excitations exhibit a mode 
softening at the BZ boundary as phases with broken translational symmetry are approached. 
We also studied the effect of NNNI on the phase diagram of the lattice bosons. We have 
proposed that such phases could be seen in dipolar condensates and discussed a number of 
ways to detect them. Recent success in achieving dipolar condensate of 52 Cr atoms and 
manipulating the interactions [1], we believe, could lead to the detection of these interesting 
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phases, particularly the SS, as well as creation of frustrated lattice bosons [23, 24] and the 
study of many associated interesting phenomena yet to be explored. 
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